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The Search for B~ 



r+z/^ at BABAR 



Luke A. Corwin (The Ohio State University) for the BaBAR Collaboration 
Stanford Linear Accelerator Center, 2575 Sand HiU Road, Menlo Park, CA 94025 

We present a search for the decay B'^ — > r^i/r using 288 fb~^of data collected at the T{4S) resonance with 
the BABAR detector at the SLAC PEP-II _B-Factory. A sample of events with one reconstructed semileptonic B 
decay {B~ l~ViX) is selected, and in the recoil a search for B'^ — > t'^Vt signal is performed. The r is 

identified in the following channels: — > e^VeVr, — » ^'^v^Vt, ■K'^Vr and — > tt^tt^Ft. We measure a 



branching fraction of B{B'* 



(0.88tQg7(stat.) ± 0.11(syst.)) x 10 * and extract an upper limit on the 



branching fraction, at the 90% confidence level, of B{B 



< 1.8 X 10 . We calculate the product of the 



B meson decay constant and |14i,| to be /s ■ = (7.0l3 6(stat.)Io g(syst.)) x 10""* GeV. 



o 



X 



1. Introduction 

1.1. Standard Model 

We here present a search for the decay 5+ 
T^Vr at the BaBAR experiment using the semilep- 
tonic tag method [T]. In the Standard Model 
(SM), purely leptonic decays of the charged B 
meson proceed via quark annihilation into a 
boson. The SM branching fraction is given by 



B{B+ -> i+iy) 



8tt 



1-^ rB\v. 



(1) 



TB, 



where Gj? is the Fermi constant, tub and m£ are 
the B meson and lepton masses, and tb is the 
B meson lifetime. The branching fraction has a 
simple dependence on Vub, the quark mixing ma- 
trix element, and /s, the meson decay constant 
which describes the overlap of the quark wave 
functions within the meson. Therefore, measur- 
ing the branching fraction could provide a clean 
experimental measurement of the SM value for 
Jb- Note that throughout this document, the 
natural units ?i = c = 1 are assumed. 

Neither this nor any other purely leptonic de- 
cay of the charged B meson has been observed. 
The most stringent currently published limit is 



B{B+ 



-Vr) < 2.6 x 10"'' at the 90% confi- 



dence level [2]. 

The SM prediction of B{B+ r+i^r) can 
be made in several different ways; we discuss 



two here. First, we substitute the experimen- 
tally determined values (when possible) and the- 
oretically calculated values (when necessary) for 
all variables into Equation [TJ All values, ex- 
cept Vub and /b, are taken from [2]. We use 
\Vub\ = (4.39±0.33) X 10-4 [3] and a lattice QCD 
calculation of /b = 0.216 ± 0.022 GeV [4]. The 
result is 

B{B+ T+Vr) = (1-59 ± 0.40) X 10""*. (2) 

The second prediction is given by the UT Fitter 
group: 

^ T+Vr) = (1.41 ± 0.33) X 10"'' [5]. 



B{B"' 



(3) 



These two values agree within uncertainty, and 
we see that current experimental limits are ap- 
proaching the SM predictions. 

1.2. Potential New Physics 

In the Type II Two-Higgs Doublet Model 
(2IIDM), purely leptonic charged B decays can 
proceed via a process identical to the SM pro- 
cess, except that the W is replaced by a charged 
Higgs boson. The relationship between the total 
and SM branching fractions is given by 



B{B+ 



B{B+ l+v) 



SM 



= 1 - tarP'P- 



•'B+ 



(4) 



where SM denotes the Standard Model branching 
fraction, tan f3 is the ratio of the vacuum expec- 
tation values of the two doublets, and rnB+ is the 
mass of the charged Higgs boson [S] . 
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2. Analysis Methods 

The data used in this analysis were coUected 
with the BaBAR detector at the PEP-II storage 
ring. The sample corresponds to an integrated 
luminosity of 288 fb"^ at the r{AS) resonance 
(on-resonance) and 27.5 fb^^ taken 40 MeV below 
BB threshold (off-resonance). The on-resonance 
sample consists of 320 x 10*^ T{AS) decays {BB 
pairs). The collider is operated with asymmetric 
beam energies, producing a boost of (3^ w 0.56 of 
the T{4:S) along the collision axis. 

The BaBAR detector is optimized for 
asymmetric-energy collisions at a center-of-mass 
(CM) energy corresponding to the T(4S') res- 
onance. The detector is described in detail in 
Ref. [7]. The components used in this analysis 
are the tracking system composed of a five-layer 
silicon vertex detector and a 40-layer drift cham- 
ber (DCH), the Cherenkov detector for charged 
TT-K discrimination, the Csl calorimeter (EMC) 
for photon and electron identification, and the 
18-layer flux return (IFR) located outside of the 
1.5T solenoidal coil and instrumented with resis- 
tive plate chambers for muon and neutral hadron 
identification. For the most recent 51 fb^^of 
data, a portion of the muon system has been 
replaced with limited streamer tubes [8j. We 
separate the treatment of the data to account 
for varying accelerator and detector conditions. 
"Runs 1-3" corresponds to the first 111.9fb~^, 
"Run 4" the following 99.7 fb"\ and "Run 5" the 
subsequent 76.8 fb~^. 

2.1. Semileptonic Tag 

In the semileptonic tagging method, one of the 
two charged B daughters of the T(4S') (here- 
after referred to as the B^) is reconstructed in a 
semileptonic decay mode B^ D'~'£~T'eX, where 
£ is e or /i. X can be nothing or a transition 
particle from the decay of a higher mass charm 
state, which we do not attempt to reconstruct. 
Our studies showed that the efficiency gained us- 
ing this method was worth the signal purity lost 
by not reconstructing the higher mass state. 

We reconstruct the candidates in four de- 
cay modes: K~tt^, K^tt^tt^tt^ , i^T^Tr+Tr", and 
K^TT^TT^ {K'^ vr+TT^). All tag i? reconstruction 



is performed assuming the only unreconstructed 
particle is the 17 g. We refer to the tracks and 
neutrals used to reconstruct the tag B as the 
"tag side" of the event. We require the net event 
charge to be zero and that all tag side tracks meet 
at a common vertex. 

2.1.1. Double Tag Sample 



Table 1 

Tag Efficiency and Systematic error derived from 
the Double Tag sample. 



Run 


^data 


Systematic 






Error 


Run 1-3 


1.05 ±0.02 


1.9% 


Run 4 


1.00 ±0.03 


3.0% 


Run 5 


0.97 ±0.03 


3.1% 



The most important control sample consists of 
"double-tagged" events, for which both of the 
B mesons are reconstructed in tagging modes, 
B- D^t-viX vs. B+ Tfi+viX. This 
sample has a very similar topology to our signal, 
and we used it to validate our Monte Carlo (MC) 
simulation and correct our tag efficiency. 

For Data and MC, we assume that the num- 
ber of reconstructed double-tagged events {N2) is 
given by N2 = e^N, where N is the total number 
of B pairs in the data sample and e is the tag 
efficiency. 

We calculated ^^^^ for Runs 1-3, 4, and 5. The 
results are shown in Table [1] The ratio is our cor- 
rection to the tag efficiency, and the fractional er- 
ror on the ratio is used as the tag side systematic 
error. 

To validate our MC, we compare double-tagged 
data events with double-tagged MC events. The 
MC samples are weighted so that the number of 
MC events generated matches the number of data 
events in the relevant sample. We find excellent 
agreement in both yield and shape between MC 
and data. 
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2.1.2. Background Rejection 

For all Z?" decay modes except X^tt+tt", the 
mass of the reconstructed D° is required to be 
within 20 MeV of the nominal mass In the 
i^T^TT+Tr" decay mode, the mass is required to be 
within 35 MeV of the nominal mass [2 . The mo- 
mentum of the tag lepton in the CM frame is 
required to be greater than 0.8 GeV. Neutral 
B mesons are excluded from the tag side by re- 
jecting events with a charged pion that could be 
combined with a to form a D*^ candidate. 

Assuming that the massless neutrino is the only 
missing particle, we calculate the cosine of the an- 
gle between the D'^£ candidate and the B meson. 



cos 9 



B~D°i 



(5) 



Here {Ejjoi^ Pd^i) a-nd (Eb, Pb) are the four- 
momenta in the CM frame, and muoi and 
are the masses of the D^t candidate and B meson, 
respectively. Eb and the magnitude of are cal- 
culated from the beam energy: Eb = i?cM/2 and 
\pb I — \/ Eg — m?g , where Eb is the B meson en- 
ergy in the CM frame. 

Correctly reconstructed candidates populate 
the range [—1,1], whereas combinatorial back- 
grounds can take unphysical values outside this 
range. We retain events in the interval —2.0 < 
cos6b-d«i < 1-1, where the upper bound takes 
into account the detector resolution and the loos- 
ened lower bound accepts those events where a 
soft transition particle from a higher mass charm 
state is missing. 

2.2. Signal Side 

After excluding the tracks and neutrals from 
the tag side, the remainder of the event (referred 
to as the "signal side" ) is searched for consistency 
with one of the above r decays, where we require 
the event contain no extra well-reconstructed 
tracks. 

Each of the four r decay modes used in this 
analysis produces exactly one charged track. The 
well-reconstructed track on the signal side is as- 
signed to a signal category based on its identi- 
fication as an e, or tt. If a signal tt can be 
combined with a tt" to yield an invariant mass 
consistent with a (see Table it is assigned 



to the T+ tt'^tt'^T't mode. 

The most important discriminating variable on 
the signal side is Ecxtra, which is the sum of all 
detected energy not associated with the tag or sig- 
nal sides of the event. For properly reconstructed 
signal events, this variable should have a peak at 
zero and a hump at higher values representing the 
X from the tag B. Thus, this analysis looks for 
excess events at low values of E'oxtra- 

We "blind" the signal region of E'oxtra in data 
until the final yield extraction is performed. In 
this analysis, the blinded region corresponds to 
-Eoxtra < 500 MeV; the region where Ecxtia > 
500 MeV is referred to as the "side band" . In 
MC, we calculate the ratio of events observed sig- 
nal region to events observed in the side band. 
We multiply this ratio by the number of data 
events observed in the side band to obtain our 
background estimate. The branching fraction and 
upper limit are extracted from the difference be- 
tween the background estimate and the observed 
number of events in the signal region. 

2.2.1. Background Rejection 

To suppress background we make requirements 
on the missing mass (Afmiss), the momentum of 
the signal candidate in the CM frame {p^ig^ad^ 
the invariant mass of an e~ and if both are 
present in our signal side (mge), the number of 
KlS in the EMC {N^j^f^i^o ) and IFR, the num- 
ber of extra tt^s in the event (iV^?*'"°), and a con- 
tinuum background suppression variable denoted 
Rt-t-. The values of this requirements are shown 
in Table [2] We do not cut on -Eoxtra; the values 
shown in the table define our signal region for 
each mode. 

Missing mass is defined as 



Mniss = \J {EriiS) - E^is)^ - (Pr(4S) - Pvis)^ 



(6) 



Here {Ey(4S)j Pt(4S)) is the four- momentum of 
the T{AS), known from the beam energies. The 
quantities -Evis and Pvis are the total visible en- 
ergy and momentum of the event which are cal- 
culated by adding the energy and momenta, re- 
spectively, of all the reconstructed charged tracks 
and photons in the event. 
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Table 2 

The selection criteria for different signal modes using a D^t^ViX tag are listed in this table. 









TT^It^Vr 


4.6 < M„„,, < 6.7 


3.2 < M™,,, < 6.1 


1.6 < 


Mrmss < 4.6 


P*signal ^ 




^■Q<P*signal 


1.7 ^ P Qig'fldl 


No IFR ifO 


No IFR K° 


No IFR if 


No IFR 


2.78 < Rrr < 4.0 


2.74 < Rrr 


2.84 < Rrr 


2.94 < Rrr 


mee > 0.1GeV/c2 








jyextra ^ 3 


jyextra < 2 


^EMCKl < 2 










selection: 
0.64 < Mp± < 0.86 GeV 
0.87 < cos 6i^_p 


i;extra < 0.31 GcV 


i;extra < 0-26 GcV 


i^cxtra < 0.48 GeV 


E,^t,, < 0.25 GeV 



Rrr is defined as 

Rrr — 



A/|"'"-0.75 GcV 
GcV 



(7) 



It is an empirically derived combination of the co- 
sine of the angle between the signal candidate mo- 
mentum and the tag i?'s thrust vector (in the CM 
frame), which is denoted cos{0rp^^ signal) ^^'^ 
minimum invariant mass constructible from any 
three tracks in an event (regardless of whether 
they are already used in a tag or signal candi- 
dates), 



which is denoted M^"". 



3. Results 

3.1. Yields 

The unblinded Data and MC are compared in 
Figure [T] and the yields for each mode are shown 
in Tabled 

3.2. Systematic Uncertainty &; Efficiency 

BaBAR has an overall systematic uncertainty of 
1.1% on the number of charged B mesons in the 
data set. A 1.5% uncertainty on the tag B yield is 
estimated using the double tag sample described 
in Section [^TTl The signal side systematic errors 
are given in Tabled) 

Overall tag B reconstruction efficiency, which 
is defined as the fraction of events tagged in a 



Table 3 

The observed number of on-resonance data events 
in the signal region are shown, together with num- 
ber of expected background events. The back- 
ground estimations include all applicable system- 
atic corrections. 



Signal T 


Expected 


Observed Events 


Decay 


Background 


in On-resonance 


Mode 


Events 


Data 


e^VeVr 


41.9 ± 5.2 


51 




35.4 ± 4.2 


36 


■K^Vr 


99.1 ± 9.1 


109 


■K^n'^T'r 


15.3 ± 3.5 


17 


All modes 


191.7 ± 11.8 


213 



BB MC sample with one B decaying to our sig- 
nal mode and the other decaying generically ac- 
cording to the branching fractions listed in [2] , is 
(6.77 ± 0.05 ± 0.10) X 10-3. In this document, 
in all cases where two uncertainties are quoted 
on a value, the first uncertainty is statistical and 
the second is systematic. The selection efficien- 
cies for each signal mode are given in Table [Sj 
where the efficiency is defined as the ratio of the 
number of signal events reconstructed (i.e. the 
rightmost column of Table [3]) to the number of 
reconstructed tagged events. 
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Table 4 

Contribution to the systematic uncertainty on the signal selection efficiencies in different selection modes. 
These uncertainties are added together in quadrature with the uncertainty on the tag B yield, extracted 
from the double-tagged control sample, of 1.5%. The uncertainty on MC statistics is added in quadrature 
to obtain the total systematic uncertainty. 



Selection 


tracking 


Particle 








Total 


Correction 


modes 


(%) 


Identification 




modeling 


modeling 


Systematic 


Factor 




0.3 


2.0 


3.6 


3.8 




5.8 


0.982 




0.3 


3.0 


3.6 


3.8 




6.2 


0.893 




0.3 


1.0 


6.2 


3.8 




7.5 


0.966 




0.3 


1.0 


3.6 


3.8 


1.8 


5.8 


0.961 



Table 5 

The signal efficiencies, mode-by-mode, relative 
to the number of tags. The branching fraction 
for the given r decay mode selected is included in 



the efficiency. 


Mode 


Efficiency (BF Included) 




0.0414 ± 0.0009 


T+ ^i+iyi_iVr 


0.0242 ± 0.0007 




0.0492 ± 0.0010 




0.0124 ± 0.0005 



3.3. Branching Fraction 

Q = ^ CL,+i,^Ps+b{Q<Qobs) (8) 

CL, = CLb = Pb{Q< Qobs) (9) 

We use a modified frequentist method, known as 
the CLs method [S], for calculating the branch- 
ing fraction and upper limit. For this method, 
we generated a large number of toy MC exper- 
iments with different branching fractions in the 
range from zero to 10 x 10"'*. 

We define the "estimator" Q, which is mono- 
tonically increasing for increasing signal, and the 
confidence levels (CL) in Equations [8] and [9l C is 
the likelihood that a given number of observed 
events could be produced only by background 
(denoted with the subscript b) or by signal and 
background (denoted with the subscript s + b). 



Pi{Q < Qobs) is the fraction of generated toy 
MC experiments that produced a Q less than or 
equal to the Q measured in data using a branch- 
ing fraction that is positive (i = s + b) or zero 
[i — b). Figure [2] illustrates our results. The left- 
hand plot shows 1 — CLfc as calculated for the 
number of events expected from different branch- 
ing fractions. Our upper limit is the branching 
fraction for which CL^ — 0.90. The right-hand 
plot shows — 21og(Q) vs. branching fraction. Our 
most likely branching fraction is at the minimum 
of this curve. The statistical uncertainty is the 
distance from the minimum where — 21og((5) is 
1.0 greater than the minimum. The resulting 
branching fraction is only inconsistent with zero 
by 1.3(7, so we present it and an upper limit, as 
well as our result for fs ■ \Vub\- 

B{B+ T+Vr) = (0.88t°:^? ± 0.11) X lO"'' (10) 
B{B+ T+Vr) < 1.8 X 10^-* @ 90% CL (11) 
fB ■ 1Kb I = {l.Otli t"oi ) X 10-4 GeV (12) 

3.4. Interpretation &; Conclusion 

Dividing the result shown in Equation ll2l bv the 
value of Vub from [3] yields Jb = O.ie^g^J^ GeV, 
which is consistent with the SM Lattice QCD pre- 
diction [4]. Our resulting branching fraction is 
consistent with both predictions shown in Sec- 
tion 11.11 (Equations [2] and [S]). Therefore, we see 
no evidence of physics beyond the SM in this anal- 
ysis. This lack of evidence can be translated into 
bounds on the properties of the charged Higgs 
boson. Figure [3] shows the regions of phase space 
that can be excluded at the 95% confidence level 
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Figure 1. i?oxtra is plotted after all cuts have been 
applied with all modes combined. MC have been 
normalized to the on-resonance luminosity. In 
addition, the background MC have been scaled 
to the data yield in the side band. Simulated 
T^Vr signal MC is plotted (lower) for com- 
parison. 



in the 2HDM. Note that mH+ < 79.3 GeV has 
been excluded by direct searches at LEP [TD]. 

A search for this mode using hadronic tags is 
underway at BABAR. We hope to publish a com- 
bined result from these two analyses in the near 
future. 
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Figure 2. The confidence level vs branching frac- 
tion is shown (left) to illustrate the extracted 
upper limit. The negative log likelihood curve 
(right) illustrates the central value and it's corre- 
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Figure 3. The shaded region shows the area ex- 
cluded by this analysis combined with the search 
by Belle for the same mode [H,. 
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